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This paper presents a further numerical study of the interaction dynamics for soh- 
tary waves in a nonhnear Dirac field with scalar self-interaction by using a fourth 
order accurate Runge-Kutta discontinuous Galerkin method. Our experiments are 
conducted on the Dirac solitary waves with a phase difference. Some interesting 
phenomena are observed: (a) full repulsion in binary and ternary collisions and 
its dependence on the distance between initial waves; (b) repulsing first, attracting 
afterwards, and then collapse in binary and ternary collisions of initially resting two- 
K^ , humped waves; (c) one-overlap interaction and two-overlap interaction in ternary 

1^ \ collisions of initially resting waves. 
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1 Introduction 



Ever since its invention in 1929 the Dirac equation has played a fundamental 
role in various areas of modern physics and mathematics, and is important for 
the description of interacting particles and fields. Consider a classical spinorial 
model with scalar self-interaction, described by the nonlinear Lagrangian L = 
iij)'-)^d^il) — mtl^il) + A(i/''j/')^ from which we may derive the nonlinear Dirac 
equation 

i^d^i^ -mil^ + 2A(V^V)V' = 0, (1) 
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where the 7^ matrices are defined by 





here a^ with /c = 1, 2, 3, denote the Pauh matrices. The nonhnear self-couphng 
term (V'V')^ i^i the Lagrangian allows the existence of finite energy, localized 
solitary waves, or extended particle-like solutions, see e.g. [1]. Several authors 
have committed themselves to analytically investigating the nonlinear Dirac 
model [2,3,4,5,6,7]. 



Here we only pay our attention to advances in numerical studies of interac- 
tion dynamics of the Dirac solitary waves. Up to now, some reliable, higher- 
order accurate numerical methods have been constructed to solve the non- 
linear Dirac equation (1). They include Crank-Nicholson type schemes [8,9], 
split-step spectral schemes [10], Legendre rational spectral methods [11], and 
Runge-Kutta discontinuous Galerkin (RKDG) methods [12], etc. The inter- 
action dynamics for the solitary wave solutions of (1) were numerically simu- 
lated in [8] by using a second-order accurate difference scheme. The authors 
saw there: charge and energy interchange except for some particular initial 
velocities of the solitary waves; inelastic interaction in binary collisions; and 
bound state production from binary collisions. Weakly inelastic interaction in 
ternary collisions is observed in [12]. The interaction dynamics in the binary 
and ternary collisions of two-humped solitary waves are first investigated in 
[13]. 

However, the experiments carried out in the literatures are all limited to the 
binary and ternary collisions of the in-phase solitary waves of (1). In this Let- 
ter we will devote ourselves to further investigating the interaction dynamics 
in the binary and ternary collisions of the Dirac solitary waves with an initial 
phase shift by using a fourth-order accurate RKDG method [12] and report 
some interesting observations. The RKDG methods adopt a discontinuous 
piecewise polynomial space for the approximate solutions and the test func- 
tions, and an explicit, high-order Runge-Kutta time discretization. It has been 
demonstrated by various experiments that the fourth-order RKDG method is 
numerically stable without generating numerical oscillation within a very long 
time interval, has uniformly numerical convergence-rates, and preserve con- 
servation of the energy and charge. We refer the reader to [12] as well as [13] 
for a detailed description and more numerical demonstrations. 



2 Preliminaries 



We restrict our attention to the (l + l)-dimensional nonlinear Dirac model (1), 
and use the same notations ps(x,t) and pQ{x,t) as ones in [12,13] to denote 
the energy and charge densities defined by 

PE{x,t) =lm{rid,^2 + r2d.^i)+m{\^,\' - l^aH - K\i'i? - \i'2?f. (2) 



pQ{x,t) =|^iP + IV'2 



(3) 



where ipi and ip2 are two components of the spinor 'ij){x,t). A standing wave 
solution of the Dirac model (1) is given as 



/ 



V'^-(x,t) 



with 



i}l'^{x,t) 



VV'r(a:,t) 




-iAt 



< A < m, 



(4) 



A{x) 



B(x) 



^j{m^ - A2)(m + A) cosh (a;^(m2 - 


-A2)) 


m + A cosh (2a;y^(m2 - A2)) 


^{{m^ - A2)(m - A) sinh [x^{m^ - 


-A2)) 



m 



A cosh {2x^{rn? - K^) 



(5) 
(6) 



The Dirac model (1) also has a single solitary wave solution placed initially 
at xq with a velocity v. 



ip''{x - Xo, t) = (V^f (x - Xo, t), V^f (x - Xo, t)) , 



(7) 



where 



^f (x-Xo,t) 
^f(x-xo,t) 



'^^r(5, i) + sign(t;) Jl-^^r (5, t), (8) 



'^V^r(5, i) + sign{v).h—^ri'"{x, i), (9) 



here 7 = 1/Vl — v"^, x = 7(x — xq — vt), i = 7(t — f (x — xq)), ^f^ and ■?/'|"' are 
defined in (4) and sign(x) is the sign function, which returns 1 if x > 0, if 
X = 0, and — 1 if x < 0. The function ■j/j''^(x — xq, t) represents a solitary wave 
travelling from left to right if f > 0, or travelling from right to left ii v < 0, 
and the standing wave i/?*"'(x — xq, t) is actually a solitary wave at rest placed 
at Xo or identical to tl)^'^{x — xq, t) with v = Q. 




Fig. 1. Dependence of pq on A and v. Left: A = 0.9; right: A = 0.1. 



The profile of the solution (4) or (7) is strongly dependent on the parame- 
ter A: it is a two-humped solitary wave with two peaks whose locations are 
determined by cos\i{2\/rnF^^ls?x) = "^^ if < A < y; it becomes a one- 
humped solitary wave with one peak located at £ = if y < A < m; and 
ijy^^{x — xo,t) = if A = m. Moreover, amplitude of the solitary waves also 
depends strongly on the velocity v: Pq{x — Xo,t) = '-ypq'^x^t). Fig. 1 shows 
that dependence, which will gives different interaction dynamics. It is worth 
noting that e^^if^^^i^x — xq, t) is still a solitary wave solution of the Dirac model 
(1), if ^ is a constant. 

In the following, our computations will work in dimensionless units, or equiv- 
alently, take m = 1 and A = |, and adopt the non- reflecting boundary condi- 
tions at two boundaries of the computational domain. The domain is covered 
by some identical cells with area of 0.05. The numerical algorithm is the P^- 
discontinuous Galerkin method in space combined with a fourth order accurate 
Runge-Kutta time discretization, see [12]. 



3 Binary collisions 



We solve (1) with the initial data 



^(a;, 0) = e''^'ilj''{x - xu 0) + e''^^ilj''{x - a:,, 0), 



(10) 



where di and Or are two real numbers, determining whether two waves are in 
phase. For convenience, we will say that two solitary waves are equal if A^ = A^ 
and \vi\ = \vr\. Throughout our numerical experiments in this section, we will 
take 9i = and 9r = n, unless stated otherwise. 




Fig. 2. The time evolution of the charge density pq. h.i = A^ = 0.5, vi 

Xr = —Xi = 10. 



0.2, 




Fig. 3. The time evolution of the charge density pq for two initially resting, equal 
waves. Left: A; = A^ = 0.6, Xr = —xi = 3, right: A; = A^ = 0.5, Xr = — x; = 6. 

3.1 Two one-humped solitary waves 



In this subsection, we study the interaction dynamics of two one-humped soh- 
tary waves with a phase shift of tt. The left plot of Fig. 2 shows the computed 
results for the case of two equal one-humped solitary waves, i.e. A; = A,, = 0.5, 
Vi = —Vr = 0.2, and Xr = —xi = 10. We see that the elastic interaction hap- 
pens and two one-humped solitary waves keep their initial shapes and velocities 
after their collisions. It is worth noting that strong overlap happens if above 
two initial waves are in phase, i.e. 6i = 6r = 0. 

When two equal waves are at rest initially and with a phase shift of tt, the 
results given in Fig. 3 show that they repulse fully each other. The repulsion 
force depends on their initial distance. The distance is smaller, the waves 
move faster in opposition. This phenomenon is different from the results on 
two waves in phase reported in [12], where the long-lived oscillating bound 
state is generated when Xr = —Xi = 3. 

Fig. 4 shows the computed results for the case of two unequal one-humped 
solitary waves. The left-hand plot is for the case of A; = 0.6, A^ = 0.8, 

—xi = 10, and the right figure is for the case 
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Fig. 4. The time evolution of the charge density pq. Left: A; = 0.6, A^ = 0.8, 
vi = —Vr = 0.2, Xr = —xi = 10; right: A; = A,. = 0.5, vi = 0.1,1;,. = —0.9, 
Xr = —xi = 10. 

of A; = A, = 0.5, vi = 0.1, fr = —0.9, and Xr = —xi = 10. We see that 
the left-hand (or right-hand) initial solitary wave transfers charge and energy 
to the right-hand (or left-hand) one and the final solitary waves are moving 
with their initial velocities. Actually, we may consider that as a traversing or 
penetration phenomenon, that is to say, two waves go through each other in 
their interaction. 

We have conducted various different experiments on binary collisions of one- 
humped waves and concluded that in general collapse phenomenon cannot be 
observed in collisions between two one-humped waves. To save space, we do 
not give corresponding plots here. 



3.2 Two two-humped solitary waves 



This subsection is to study the interaction of two two-humped solitary waves 
with a phase shift of tt. 

Fig. 5 shows the computed results for the cases of A; = A, = 0.1, f; = —Vr = 
0.2 and 0.9, and Xr = —xi = 10. We observe that the final solitary waves keep 
their initial velocities but with different shapes; the collapse does not happen 
in these both cases. It is worth noting that the collapse will happen if those 
two waves are in phase, see Fig. 2 in [13]. 

Consider the case of that two unequal waves are at rest initially and with a 
phase shift of tt. The results given in Fig. 6 show that they repulse each other 
essentially, and the repulsion force depends on their initial distance. When two 
initial waves stand more nearly, the repulsion dominates in their interaction, 
thus they move outside fast and cannot re-collide each other, see the left figure 
in Fig. 6. But when the distance is relatively big, the right plot of Fig. 6 shows 
that two waves first repulse each other, then attract afterwards and collapse. 
When two initial waves are equal, at rest, and with a phase shift of n, we only 




Fig. 5. The time evolution of the charge density pq. A/ = A^ = 0.1, Xr 
Left: vi = —Vr = 0.2; right: vi = —Vr = 0.9. 



-XI = 10. 




Fig. 6. The time evolution of the charge density pQ for two initially resting waves. 
A^ = 0.1, Ar = 0.2. Left: Xr = —xi = 4; right: Xr = —xi = 5. 

observed full repulsion which is similar to that in the left figure of Fig. 6. 

Besides the above, we also conduct some other experiments and find that 
collapse happens easily in collisions of two in-phase, equal, two-humped waves, 
but it may not appear in collisions between two in-phase, unequal, two-humped 
waves, or two equal two-humped waves with a phase shift of n, or two unequal, 
two-humped waves with a phase shift of n. 



3.3 A one-humped solitary wave and a two-humped solitary wave 



This subsection is to study the interaction of a one-humped solitary wave and 
a two-humped solitary wave, which are with a phase shift of n. 

Fig. 7 shows the computed results for the case of A; = 0.1, A^ = 0.9, vi = 
—Vr = 0, and Xr = —xi = 6. We see the quasi-stable long-lived oscillating 
bound state, which is essentially same as one shown in Fig. 3 in the paper 
[13]. Generally, when there is a big difference between the peak values of two 
initial waves, macroscopical behavior of the interaction dynamics are essen- 
tially independent on their phase shift. 




Fig. 7. Bound state formed in binary collisions between two initially resting, 
out-of-phase waves. A; = 0.1, A^ = 0.9, vi = Vr = 0, and Xr = —xi = 6. 




Fig. 8. The time evolution of the charge density pq. A/ 
Xr = -xi = 10. Left: A^ = 0.5; right: A^ = 0.9. 
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0.2, 



Fig. 8 shows the computed results for the cases of A^ = 0.1, A^. = 0.5 and 0.9, 
vi = —Vr = 0.2, and Xr = —xi = 10. It tells us that collapse may be observed 
in collisions between a one-humped solitary wave and a two-humped solitary 
wave when they do initially travel face to face, but it may also not happen. 



4 Ternary collisions 



In this section, we study ternary collisions by solving the nonlinear Dirac 
model (1) with the following initial data 



je^ 



\9r „I.SS I 



V(a;, 0) = e'^il^Ti.^ - ^u 0) + e''-V::(a: - x^, 0) + e''^il^l\x - x,, 0), (11) 









Fig. 9. Collisions of three initially resting two-humped waves: A/ = A^ = A,, = 0.1, 
Xr = —xi = 10, Xm = 0, 01 = Or = 0, 6m = VT. Left: plot of pQ{x,t); right: plots of 
PQ{0,t) and pEiO,t). 



where 6i, dm and Or are three real numbers, determining the initial phase of 
corresponding waves. 

The first case we consider is collisions of three initially resting two-humped 
solitary waves: A/ = Km = A^ = 0.1, Vi = Vm = v.^ = 0, Xr = —Xi = 10, 
Xm = 0, 9i = 9r = 0, and 9m = tt. The results are shown in Fig. 9. We see from 
the left plot that the waves first repulse each other because two neighboring 
waves are out-of-phase, but after t = 100 they begin to attract each other, and 
final collision results in collapse. Symmetry of the solutions is kept very well. 
The right figure of Fig. 9 gives the charge and energy densities at a; = as a 
function of time. We observe that before collapse happens, the middle wave is 
oscillating because of bind from left and right waves although its displacement 
seems to be unchanged. 



The second case is collisions of three initially resting one-humped waves: A; = 

Am = Ar = 0.5, Vl = Vm = Vr = 0, Xr = —Xl = 10, Xm = 0, 6'; = 71, and 

9m = 9r = 0. The results are shown in the left plot of Fig. 10. We see that 
the first interaction happens between two right in-phase waves around t = 45, 
and then two faster moving waves are formed. The initially resting left wave 
begin to be moving towards left due to repulsion between it and the right 
waves, and then it is catched up with and interacted by the left moving wave 
generated in the first interaction around t = 130. Overlapping happens in all 
two interactions. If we consider collisions of three initially resting one-humped 
waves with A; = A^ = A^ = 0.6 or 0.9, the interaction does only happen 
between two right neighboring in-phase waves, see the right plot of Fig. 10. 
The reason is that the left-moving (middle) wave formed in the interaction of 
two initially in-phase (right) waves is not faster than the left wave. 




Fig. 10. The time evolution of the charge density pq. vi 
Xr = —xi = 10, Xm = 0, Oi = vr, 9m = Or = 0. Left: A^ 



A„ 



K 



-- 0, 
0.5; 



right: A; = A^ = A^ 



0.6. 



5 Discussions and Conclusions 



In this paper, we further studied the interaction dynamics for the solitary 
waves of a nonlinear Dirac model (1) with scalar self-interaction. By using a 
fourth-order accurate RKDG method presented in [12], we investigated the 
interaction of the Dirac solitary waves with an initial phase shift of tt for 
the first time, and observed that: (a) full repulsion in binary and ternary 
collisions, and the initial distance between waves is smaller the repulsion is 
stronger; (b) repulsing first, attracting afterwards, and then collapse in binary 
and ternary collisions of two-humped waves; (c) interaction with one overlap 
and two overlaps in ternary collisions of initially resting waves, which depends 
on initial parameter A. We concluded that in general collapse phenomenon 
cannot be observed in collisions between one-humped waves, but it happens 
easily in collisions of in-phase, equal, two-humped waves; the macroscopical 
behavior of the interaction dynamics is essentially independent on their initial 
phase difference, when there is a big difference between the peak values of 
initial waves. Although we have investigated the influence of initial phase 
difference on the interaction dynamics of the Dirac solitary waves, it will be 
interesting and important to study the evolution of the relative phase of the 
solitons during their collisions. 
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